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We evaluate the contribution of a class of higher-twist operators to the lowest moment of the Structure Func- 
tions, by computing appropriate matrix elements of six four-fermion operators in the quenched approximation. 
Their perturbative renormalization constants and mixing coefficients are calculated in the 't Hooft-Veltman 
scheme of dimensional regularization, using codes written in the algebraic manipulation computer language 
FORM. 



1. INTRODUCTION 

By performing analytic calculations and 
numerical simulations with quenched Wilson 
fcrmions, we have computed the contributions 
of some classes of higher-twist operators to mo- 
ments of the structure functions. 

These moments are related via an OPE to 
hadron matrix elements of local operators, whose 
twist is the difference between their dimensions 
and their spin. Twist-2 operators give the lead- 
ing contributions, while twist-4 operators corre- 
spond to the power corrections: 



A^n(g') = J dxx"-2_p2(x,Q2 

= c^'\Q')^,\g{^l))A^\^,) 
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The leading twist contribution can be written as 

= E/ Q) (x/"^): where Q/ is the charge of 
the quark of flavor /. At twist-4 level we consider 



*Talk given by S. Capitani at Lattice 99, Pisa (Italy). 



the 1 = 2 pion structure function ||l| 

^ F^^ + Ff - 2Ff , (2) 

which belongs to a flavor 27-plet, receives contri- 
butions only from 4-fermion operators and there- 
fore cannot mix with twist-2 operators □. While 
the Wilson coefficients C^"^ can be calculated in 
perturbative QCD the computation of the 

matrix elements is a strong interaction prob- 
lem that cannot be treated perturbatively. Some 
results have been obtained in the past using the 
MIT bag model , but a calculation in a reliable 
model-independent way and from first principles 
is only possible in lattice QCD. 

The leading-twist contribution to the lowest 
unpolarized moment is given by the operator 

= ^V'Tmi D^^2 i' ~ traces, (3) 

^When mixing with operators of leading twist is forbidden, 
one gets rid of renormalon ambiguities completely. When 
such a mixing occurs, a non-perturbative comnutation of 
the Wilson coefficients in the way proposed in [EIH] would 
avoid renormalon problems. That method will be a useful 
complement to the present one for an understanding of 
higher— twist effects. So far however only 2— quark states 
have been simulated. 



2 



with Wilson coefficient cf^ = 1 + 0(g2) 
and forward matrix elements (pIO^^j^^^} If*) ^ 

(2) 

2A2 \PfiiPfj.2 ~ traces]. TheLtwist-4 contribution 
involves various operators a, but only 4-fermion 
operators can contribute to the 1 = 2 combina- 
tion (|^), and the only 4-fermion operator that 
appears in at order g-^ is 

Alu^2 = ^It^il^t"^ 'iplt.2l5t""ip - traces, (4) 

with Wilson coefficient C^'^^ = .92(1 + 0(32)) 
and forward matrix elements (pI^i^j^,}. |p) = 

2A^^ [Pfj.iPp2 —traces] . This operator mixes under 
renormalization, and on the lattice there is even 
more freedom to mix, given the lower symmetry 
of the theory compared to the continuum. 

2. PERTURB ATI VE 

RENORMALIZATION 

We have computed pion matrix elements of 
six 4-fermion operators and perturbatively calcu- 
lated their renormalization factors, which relate 
the lattice numbers to physical quantities in the 
MS scheme, using the 't Hooft-Veltman prescrip- 
tion for 75. We consider the following operators 
in Euclidean space, symmetrized in fi and ly; 





= tpjfj.t"'^ ^ji.t"'^ — traces 


(5) 




= i'lti.lbt"'^' V'Ti^Ts^"^ - traces 


(6) 


rpC 


— 4"^p.p't°'4' 4"^pi't°'''P ~ traces 


(7) 




= ipjpip ip^i^ip — traces 


(8) 




= i^lplzip i^lvlbi^ - traces 


(9) 


T 

pu 


— V'O'mp''/' '4"^pu'4' ~ traces. 


(10) 



Not all above operators are present in the 
OPE (|^), they appear however wheruone com- 
putes the lattice radiative corrections 

One needs the renormalization factors of the 
operators on the lattice as well as in the contin- 
uum, in order to obtain physical continuum ma- 
trix elements from the Monte Carlo simulations. 

^For discussions of twist-4 operator bases and general 
higher-twist effects, see also |pU2lp|{. 

*In principle also gauge— variant operators could appear in 
the mixing, but there are no such 4— fermion operators with 
dimension 6, and 2-fermion operators do not contribute. 



The connection between the two is 
where 

j 

are the renormalized (in our case in the MS 
scheme) continuum and unrenormalized lattice 
1-loop expressions. We computed the matrix 
elements on quark states, in a general covariant 
gauge. The differences AR,j = R\f - i?™"* enter 
then in the renormalization factors 

Zij (afi, g) = 5^J " -^^l^Rij , (13) 

which convert from the lattice to the continuum, 
and are the object of our computations. 

To evaluate the Feynman diagrams and ob- 
tain the algebraic expressions for the renormal- 
ization factors, we run computer codes written in 
the symbolic manipulation language FORM. To 
properly deal with the 75 matrices we have also 
written additional computer routines which are 
able to perform computations with the 't Hooft- 
Veltman prescription. This is the only one proven 
to be consistent for 75, and if we were not using 
it we would not find the right mixing factors be- 
tween axial, vector and tensor operators. 

Implementing the 't Hooft-Veltman scheme in 
a computer program is a challenging task, com- 
pounded by the fact that we are dealing here 
with lattice calculations. In fact, d-dimensional 
sums are split in 4 and {d — 4) dimensions, in 
which the Dirac algebra is different from usual 
(see for example [Q). The treatment of 7- 
matrices, already complicated by the non-validity 
of the summation convention on the lattice (as 
Lorentz invariance is broken), is then subject to 
much more complicated rules. When running 
the codes, the sum splitting increases the num- 
ber of terms by about one order of magnitude, 
and this slows down considerably the computa- 
tions, besides generating problems in memory 
management as well. We have found the compu- 
tations presented here, which only concern Wil- 
son fermions without any improvement, to be al- 
ready very demanding. 
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We use the Kawai method Q in deaUng with di- 
vergences, which allows a strong check in the ana- 
lytic lattice calculations by using two different in- 
termediate regularizations: dimensional regular- 
ization with the 't Hooft-Veltman prescription on 
the lattice, and an infrared-mass regularization. 

Using the 't Hooft-Veltman prescription for the 
75 matrices, the one-loop renormalized operators 
in the MS scheme can be written ip terms of the 
lattice operators (for /i = 1 /a) as 



tion to the lowest moment of 



c\R 



V - (^0.281578 V + 0.000532 A 
-h0.000997 A'- + 0.022899 T 

V -g^ (o.348170 V + 0.002393 A 

A" ~ 5^1^0.291756 A" + 0.000532 V 
-H0.000997 V + 0.006785 T 
-h0.012722 T'^) 

A - .g2 (^0.266750 A + 0.002393 V 
+0.030533 T'^) . (14) 



The complete results for general Nc are given 
in px| ]. Although used here for the pion, they 
are of more general application, and we will use 
them also for the nucleon case. 



3. RESULTS 

We made simulations on a 16^ X 32 lattice 
at quenched (3 = 6.0, on a Quadrics QH2 at 
DESY Zeuthen, with a statistics of 400 config- 
urations. To extract the expectation values of 
the operators we computed ratios of three-point 
to two-point functions of the pion, Ro{t,T) — 
{n{t)O{T)ni0))/{n{t)n{0)). 
After extrapolating the matrix elements for 
to the chiral limit and using the renormal- 
ization factors (O) at = /i^ = a^^, wc obtain 



A^'') = 0.085(27)/2, which multiplied by C^'^^ and 
the kinematical factor gives the twist-4 contribu- 



^For a similar calculation in the context of weak matrix 
elements, involving a different set of operators and a dif- 
ferent scheme, see also ||. 



M. 



(4) 



0.085(27) 



2„2 



1^9 



(15) 



For the leading contribution we have (x^) = 
0.273(12) (independent of flavor), so that weight- 
ing with the charges, say of the 7r+, one gets 11 

(16) 



Mf^""^ = 0.152(7) -fO(g2). 



The relative magnitude of the twist-4 and twist- 
2 contributions to the lowest moment of the pion 
structure functions turns out to be small, but is 
expected to be larger for the nucleon, where the 
scale /tt is replaced by . We are now studying 
the nucleon with the same machinery and using 
the same renormalization factors as here. 
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